Plastic flow, voltage bursts, and vortex avalanches in superconductors 
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We use large-scale parallel simulations to compute the motion of superconducting magnetic vortices 
during avalanches triggered by small field increases. We find that experimentally observable voltage 
bursts correspond to pulsing vortex movement along branched channels or winding chains, and relate 
vortex flow images to features of statistical distributions. As pin density is increased, a crossover 
occurs from interstitial motion in narrow easy-flow winding channels with typical avalanche sizes, to 
pin-to-pin motion in broad channels, characterized by a very broad distribution of sizes. Our results 
are consistent with recent experiments. 
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Recent studies of systems characterized by avalanche 
dynamics, where energy dissipation occurs in sudden 
bursts of collective activity, have revealed many interest- 
ing dynamical behaviors. Granular assemblies Q, mag- 
netic domains [0, fluid flow 0, and flux lines in type-II 
superconductors [||-|| all display avalanches, but produce 
broad distributions of avalanche sizes only under cer- 
tain conditions, indicating a lack of universal response to 
small perturbations. Probing the microscopic origins of 
different avalanche behaviors by experimentally record- 
ing the microscopic motion over time scales long enough 
to statistically characterize the system is difficult. In 
contrast, numerical simulations provide both precise con- 
trol of microscopic parameters and dynamical informa- 
tion, making them a useful tool. Early avalanche sim- 
ulations used simple discrete models [61, but recent ad- 
vances in parallel processing allow the use of more realis- 
tic continuous molecular dynamics (MD) models. In this 
paper, we present extensive parallel MD simulations of 
flux-gradient-driven superconducting vortex avalanches. 
By varying the pinning density n p and maximum pin- 
ning strength f p , we obtain avalanche distributions in 
good agreement with recent experiments JiJjTj , and quan- 
tify how the density and strength of pinning sites af- 
fects the breadth of these distributions. Characteristic 
avalanche sizes and lifetimes appear at low pin densities 
when narrow easy-flow winding channels of interstitial 
vortices form. At higher pin densities, pin-to-pin trans- 
port through vortex chains occurs, and the avalanche size 
and lifetime distributions remain very broad. We find 
that there is no universal distribution valid for all pin- 
ning parameters. 

Flux penetrates a type-II superconductor in the form 
of discrete quantized vortices. As an external field is 
slowly increased, a metastable gradient in vortex density, 
termed the Bean state Jlo[ , forms as vortices are driven 
into the sample by their mutual repulsion and are held 
back by defects in the material . To model this sys- 
tem, we simulate an infinite slab with a magnetic field 



H = Hi applied parallel to the surface so that there 
are no demagnetization effects. The rigid vortices and 
straight columnar pins we consider are all parallel to z, so 
we can obtain all relevant dynamical information from a 
transverse two-dimensional slice, in the x-y plane, of the 
three-dimensional slab. The flux lines evolve according 
to a T = MD algorithm. The vortex-vortex repulsion, 
given by the modified Bessel function Ki(r/X), is cut off 
beyond r = 6A, where A is the penetration depth, so that 
each vortex interacts with up to 100 neighbors, and im- 
portant collective effects, neglected in simulations with 
shorter interaction ranges, are observed. Each 24A x 26 A 
sample contains up to 3700 attractive parabolic pins of 
radius £ p = 0.15A with pinning densities n p — 0.96/A 2 , 
n p = 2.40/A 2 , or n p = 5.93/A 2 , and pinning strengths 
uniformly distributed over the range /™ ax /5 to /™ ax , 
where /™ ax = 0.3/o, l-0/o, or 3.0/o- All forces are given 
in units of /o = <I>o/87r 2 A 3 and lengths in units of A. 

The total force on vortex i is given by fi = f ™ + f" p = 
rjvi, where the force on vortex i from other vortices is 
f™ = Z)i=i fo K i(\ T i - Tj\/X)rij, and the force from 



pinning sites is = E£i(/p/£p)l r * ~ r k'\ 9 ((&> ~ 
|rj — rj^|)/A)?jfc. Here, 8 is the Heaviside step func- 
tion, Yi (vi) is the location (velocity) of vortex i, r!l 
is the location of pinning site k, there are N p pin- 
ning sites and N v vortices, fy = (rj — r 3 -)/|rj — r,-|, 



» 



)/h 



and we take r\ = 1. Us- 



ing roughly 10 4 hours on an IBM SP parallel computer, 
we recorded more than 10 4 avalanches for each of five 
combinations of n p and f p . In addition to the results de- 
scribed in this paper, we obtained distributions of vortex 
velocities and displacements as well as analytical results 
supporting our observations [Tl| ] . Further simulation de- 
tails appear in 

A slowly increasing external field is modeled by adding 
a single vortex to an unpinned region along the sample 
edge whenever the system reaches mechanical equilib- 
rium [jll] [l3| . The majority of these small field increases 
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FIG. 1. Consecutive snapshots of the velocity field of the 
vortex lattice in a 25A x 5A region of a 26A x 24A sample 
with a high density of strong pinning sites: n p = 5.93/A 2 , 
ymax _ 3 qJq. Each moving vortex is indicated by an arrow 
with its length scaled by the velocity of the vortex; the re- 
maining vortices are indicated by small crosses. All vortices 
sit in pinning sites (not shown here) when not in motion. The 
disturbance propagates from the dense left edge of the sam- 
ple to the relatively less dense right edge. The vortices in the 
rest of the sample (not shown) were not depinned. The time 
interval between snapshots is the typical time th required for 
a vortex to hop from one pinning site to another. The illus- 
trated motion is typical for a medium-size avalanche in this 
sample. 

result in only slight shifts in vortex positions, but oc- 
casionally one or more vortices will become depinned, 
producing an avalanche. We find that avalanche distur- 
bances propagate as an uneven pulse, as seen in Fig. [j] 
where arrows indicate instantaneous vortex velocities. 
Events with longer lifetimes often contain more than one 
pulse of motion (i.e., multiple oscillations in the total 
avalanche velocity) (jf3. 

By imaging individual avalanches in our samples, we 
find that a chain of vortices is displaced in a typical 
event. Each vortex in the chain is depinned, moves a 
short distance, and comes to rest in a nearby pinning 
site. Vortices outside the chain transmit stress by shifting 
very slightly inside pinning sites, but are not depinned. 
In Fig. ^, filled circles represent the initial positions of 
vortices that were depinned during the avalanche, while 
small crosses mark the vortices that remained pinned. 
Easily visible chains of moving vortices extend down the 
flux gradient in the x-direction, winding slightly in the 
y-direction. Chain size varies from event to event: in 
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FIG. 2. Snapshots of vortices depinned in an avalanche for 
two different events in the same sample. Here, stationary 
vortices are represented by crosses, while the initial positions 
of vortices that were depinned are marked with filled circles; 
the pinning sites are not shown. Vortices move towards the 
right of each panel, in the direction of decreasing vortex den- 
sity. The entire 26A x 24A sample, with a pin density of 
n p = 2.40/A 2 and /™ ax = 3.0/ , is shown. 

some cases a chain spans the sample, as shown in Fig. |^, 
while in other events the chain contains only three or 
four vortices. In each case, although the disturbance may 
cross the sample, an individual vortex does not. Thus, 
the time span of a typical avalanche is much shorter than 
the time a single vortex takes to traverse the sample. 

When the pinning density is high, chains of moving 
vortices are equally likely to form anywhere in the sam- 
ple. As the pinning density is lowered, vortices move only 
through narrow well-defined interstitial winding channels 
in which vortices are weakly held in place only by the re- 
pulsion of other vortices which sit in pinning sites. To 
identify the cumulative pattern of flow channels for dif- 
ferent pinning parameters, in Fig. ^ we plot vortex tra- 
jectories with lines over the course of many avalanches. 
A concentration of trajectory lines indicates a heavily- 
travelled region of the sample. In Fig. [|(a) we show 
a small portion of a sample with low pinning density. 
We see that all motion occurs through narrow easy-flow 
interstitial channels. In these channels, mobile intersti- 
tial vortices move plastically |l5| around their strongly 
pinned neighbors in a manner similar to that recently 
imaged experimentally ||. This sample also contains 
vortices oscillating in interstitial "magnetic traps." As 
the pin density increases, the amount of interstitial pin- 
ning decreases and the number of channels increases, as 
in Fig. ||(b), until at high pin densities [Fig. ||(c)], there is 
no interstitial flow and the vortices move only from pin to 
pin. Here, where no easy-flow channel exists, avalanches 
are spread evenly throughout the sample. 

We determine how the microscopic pinning parameters 
affect avalanche size by finding the total avalanche life- 
time r for each event. Since vortices typically hop from 
pin to pin in samples with high pinning density, a nat- 
ural unit of time is the interval th a vortex spends hop- 
ping between pinning sites, and so we use scaled lifetimes 
t* = r/th- To find th, we assume that 
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FIG. 3. Lines indicate the paths vortices follow through a 

portion of the sample over an extended period of time cover- 
ing many avalanches. Here, filled dots represent vortices, and 
open circles mark pinning sites. Each panel shows an 8A x 8A 
region of 24A x 26A samples with /™ ax = 3.0/o and differ- 
ent pinning densities: n p — 0.96/A 2 (a), n p — 2.40/A 2 (b), 
and n p = 5.93/ A 2 (c). The presence or absence of easy-flow 
channels is clearly dependent on pin density. The channels 
present in (a) and (b) lead to avalanches with characteris- 
tic sizes and lifetimes superimposed on a broad distribution. 
Samples with higher pinning density produce very broad dis- 
tributions of avalanche sizes. 



each vortex hops a distance d p 



-1/2 



tance between pinning sites, and that the vortex speed 
v c is proportional to the depinning force, v c — |fj|/?7, 
where |^| w -f p . This gives t h = d p /v c ps rjf^rip 1 ' 2 . In 
Pbln, for example, 77 ~ 0.3 G-m/f2, so using f p — 3.0/o, 
n p = 5.93/A 2 , and A ~ 65 nm gives th ~ 150 ns. The 
plot of P(t*) in Fig. |](a) for samples with dense pinning 
shows that in each case the distribution is very broad 
and can be written as P{t*) ~ ( T *) -1 ' 4 over a range 
of r*. If the pinning density n p is reduced, the form 
of P{t*) changes noticeably, as in Fig. ||(b). Here, we 
find an enhanced probability for avalanches with a char- 
acteristic value as a result of the appearance of easy- 
flow interstitial channels. Many avalanches in these sam- 
ples consist of a single sample-spanning pulse of motion 
through one of these channels. The estimated pulse life- 
time produced by a straight channel is r cst ~ thNh, 
where Nh is the number of vortices in the channel. Since 
Nh ~ L x ^/n^, where L x is the sample length, we find 
Test ~ vLxfp 1 \Jn v /n p , which gives a characteristic value 

Test/th « L X V^ = 26AV1-5/A 2 « 30. This 

agrees well 



of t* = T cst /t h w L xy ^ = 26Aa/1.5/A 2 
value, indicated by an arrow in Fig. |J(b), 
with the peak in the distribution of r*. 

For all pin densities, only a small fraction of the vor- 
tices move significantly while the rest shift in pinning 
sites, as seen by considering the distance di each vor- 
tex is displaced in an avalanche. Vortices that hop from 
pin to pin create a peak in P(rfj), marked with arrows 



in Fig. |4|(c-d), at di « d p — rip 1 ^ 2 . For those vortices 
that remain pinned and accumulate stress in the vortex 
lattice, di < d p , we can approximate the distribution 
by P(di) ~ d^ p , where p ~ 1.4 for all samples except 
p ~ 1.2 for [n p = 0.96/A 2 , /™ ax = 3.0/ ], and p ~ 0.9 for 
[nv = 5.93/A 2 , /™ ax = 0.3/o]. An analytical argument 
jyj, sketched here, predicts a similar p for all samples 
since the distribution is generated only by pinned 
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FIG. 4. Avalanche distributions: lifetimes r* (a-b); in- 
dividual vortex displacements di (c-d); total number 7V a of 
vortices moving during an event (e-f); number Nf of vor- 
tices falling off the edge of the sample (g-h). The left panels 
(a,c,e,g) refer to samples with high pin density n p — 5.93/A 2 , 
and differing pinning strengths: filled triangles (solid line), 
/™ ax = 3.0/o; open diamonds (dashed line), /™ ax = l.O/o; 
plus signs (heavy solid line), /™ ax = 0.3/o. The right pan- 
els (b,d,f,h) refer to samples with /™ ax = 3.0/o and varying 
pinning densities: filled triangles (solid line), n p = 5.93/A ; 
open squares (dot-dashed line), n p = 2.40/A 2 ; asterisks (dot- 
ted line), n p = 0.96/A 2 . In (a) the hopping time is th = 6, 
18, and 57 MD steps, respectively, while in (b), th — 6, 17, 
and 53 MD steps, respectively. 

vortices and is not affected by the presence of easy-flow 
channels. The addition of a vortex to the sample ex- 
erts a small additional force / on an arbitrary vortex 
located a small distance r away (r <C A), displacing this 
vortex a distance Su(r) = {f/rj)8t ~ K\{r/\) ~ 1/r. 
Since there are 5N(r) = 2-nm v br vortices a distance r 
from the added vortex, we find p = —dhi6N/dhx5u — 
— (l/r)(— l/r) _1 = 1, in general agreement with our com- 
puted values of p, p « 0.9 - 1.4. 
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Pinning density n p has a significant impact on the 
number of vortices N a that move a distance di larger 
than the pin diameter 2£ p . We plot P{N a ) in Fig. ||(e- 



f). When n p is low, 



< 



2.4/A , vortices move only in 



narrow easy- flow winding channels, producing a charac- 
teristic value of N a , N a fa Nh, marked with an arrow in 
Fig. |i|(f). The strongly pinned vortices around the in- 
terstitial channel cannot be depinned, so large values of 
N a are not observed. At high pinning densities, as in 
Fig. ^(e) where n p — 5. 93/ A 2 , we can approximate the 
form of P(N a ) for small N a as P(N a ) ~ N~P. As the 
pin strength decreases, decreases slightly: ~ 1.4 for 
/ p max = 3.0/d, ~ 1.0 for = l.O/o, and ~ 0.9 

for /™ ax = 0.3/o- This is in contrast to the lifetime 
distributions which were not affected by changing pin- 
ning strength. Smaller values of indicate a relative 
increase in the frequency of large avalanches compared 
to small ones. The trend in indicates the impor- 
tance of avalanche width in determining N a . In sam- 
ples with strong pinning, avalanche width is suppressed 
since strongly pinned vortices on either side of the mov- 
ing chain are not depinned. As the pinning weakens, 
wider avalanches occur when weakly pinned vortices ad- 
jacent to a slowly moving chain depin and join the mo- 
tion. Thus, avalanches with similar scaled lifetimes are 
likely to be wider in samples with weak pinning than in 
samples with strong pinning. 

Altering the pinning parameters affects the number of 
vortices Nf that exit the sample during an avalanche, as 
shown in the plot of P(N f ) in Fig. |(g-h). If we ap- 
proximate P(Nf) for low N f by the form P(N f ) ~ Nj a , 
we find that all samples with high pinning density, n p = 
5.93/A 2 , have a ~ 2.4. As n v decreases, a increases: 
a ~ 3.4 for n p = 2.40/A 2 and a ~ 4.4 for n p = 0.96/A 2 . 
When all vortex motion occurs in an easy-flow intersti- 
tial channel, a increases since the channel does not build 
up enough stress to allow events with large Nf to oc- 
cur. Smaller events continually relieve the accumulated 
stress instead. For example, with the small number of 



flux paths in Fig 
and a ~ 3.4. Fig. 



(b), events with large Nf are rare, 
(a) corresponds to the extreme case of 
a sample with only one channel, for which n p = 0.96/A 2 
and a ~ 4.4. In samples with high pin density, even after 
the stress in one vortex path has been depleted by a large 
avalanche, other areas still contain enough stress to re- 
main active in large and small events while the depleted 
regions build up stress again. This leads to a greater like- 
lihood of large events and correspondingly smaller values 
of a, a < 2.4. Distributions similar to this high pin den- 
sity case have been obtained experimentally in pj , where 
values of a ranging from 1.4 to 2.2 are observed. The pin- 
ning density from grain boundaries in the experimental 
sample is very high, n p ~ 100/ A 2 , so it is reasonable that 
the a values observed in Q] are similar to the a values 
produced by our most densely pinned samples. Broad 
distributions with a's of 1.7 to 2.2 were also observed in 



0, in good agreement with both Q and our results. In 
addition, (t]] finds a regime where avalanches of a char- 
acteristic size occur, offering an experimental example in 
which samples with a lower density of weaker pinning 
sites produce narrower distributions, as we also observe. 

We have quantitatively shown how pinning determines 
the nature of vortex avalanches. By using large-scale 
MD simulations, we image pulses of motion in chain- like 
disturbances through the sample. The presence or ab- 
sence of distinct channels for flow leads to a crossover 
from broad distributions of avalanche size to character- 
istic sizes. Pinning strength causes a transition between 
strongly plastic flow to mildly plastic, "semi-elastic" flow. 
Lowering the pinning density causes a transition from 
broad distributions to distributions with characteristic 
sizes. 
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